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1 Introduction 



Let N be an integer equal or greater than 2 and k„ = (_i)i+^/2 jf ]\[ jg even, k^^ = ±1 if N is odd. 
Consider the heat-type equation of order N: 

du d^u 

For N — 2, this equation is the classical normalized heat equation and its relationship with linear 
Brownian motion is of the most well-known. For > 2, it is known that no ordinary stochastic process 
can be associated with this equation. Nevertheless a Markov "pseudo-process" can be constructed by 
imitating the case N — 2. This pseudo-process, X — (X{t))t>o say, is driven by a signed measure as 
follows. Let p{t;x) denote the elementary solution of Eq. p.ip . that is, p solves p.ip with the initial 
condition p{0; x) — 6{x). This solution is characterized by its Fourier transform (see, e.g., [13]) 

+00 
'OO 

The function p is real, not always positive and its total mass is equal to one: 

p{t; x) dx = 1. 



Moreover, its total absolute value mass p exceeds one: 

/+00 
\pit;x)\dx > 1. 
-oo 

In fact, if N is even, p is symmetric and p < +oo, and if N is odd, p — +oo. The signed function p is 
interpreted as the pseudo-probability for X to lie at a certain location at a certain time. More precisely, 
for any time t > and any locations x,y € M, one defines 

F{X{t) e dy\X{0) - x}/dy = p{t; x - y). 

Roughly speaking, the distribution of the pseudo-process X is defined through its finite-dimensional 
distributions according to the Markov rule: for any n > 1, any times ti, . . . ,tn such that < < • • • < i„ 
and any locations x, t/i, ...,?/„ G M, 

n 

P{X{h) e dyi,...X(i„) e dyn\X{0) ^ x}/dyi . . .dyn =Y[p{U - U_i;y,_i - y,) 

i=l 

where = and yo = x. 

This pseudo-process has been studied by several authors: see the references [2] to [4] and the references 
[8] to [20]. 

Now, we consider the sojourn time of X in the interval [0, +oo) up to a fixed time t: 

T{t)= f l[o,+oo)(^(s))ds. 
Jo 

The computation of the pseudo-distribution of T{t) has been done by Beghin, Hochberg, Nikitin, Ors- 
ingher and Ragozina in some particular cases (see [2j H] E] [HI [20]) , and by Krylov and the second author 
in more general cases (see [TUl fTT] ). 

The method adopted therein is the use of the Feynman-Kac functional which leads to certain dif- 
ferential equations. We point out that the pseudo-distribution of T{t) is actually a genuine probability 
distribution and in the case where N is even, T{t) obeys the famous Paul Levy's arcsine law, that is 

P{T(t)eds}/ds= 



■7Ty/s{t - s) 



We also mention that the sojourn time of X in a small interval (— e, e) is used in [3] to define a local time 
for X at 0. The evaluation of the pseudo-distribution of the sojourn time T{t) together with the up-to- 
date value of the pseudo-process, X{t), has been tackled only in the particular cases = 3 and iV = 4 
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by Beghin, Hochberg, Orsingher and Ragozina (see OS]). Their results have been obtained by solving 
certain differential equations leading to some linear systems. In [3J |31 [TT], the Laplace transform of the 
sojourn time serves as an intermediate tool for computing the distribution of the up-to-date maximum 



In this paper, our aim is to derive the joint pseudo-distribution of the couple (T{t),X{t)) for any 
integer N. Since the Feynman-Kac approach used in (21 14] leads to very cumbersome calculations, we 
employ an alternative method based on Spitzer's identity. Since the pseudo-process X is properly defined 
only in the case where N is an even integer, the results we obtain are valid in this case. Throughout the 
paper, we shall then assume that N is even. Nevertheless, we formally perform all computations also in 
the case where TV is odd, even if they are not justified. 

The paper is organized as follows. 

• In Scction[2l we write down the settings that will be used. Actually, the pseudo-process X is not well 
defined on the whole half-line [0,-|-oo). It is properly defined on dyadic times A;/2", k,n gN. So, 
we introduce ad- hoc definitions for X{t) and T{t) as well as for some related pseudo-expectations. 
For instance, we shall give a meaning to the quantity 



which is interpreted as the 3-parameters Laplace-Fourier transform of {T{t),X{t)). We also recall 
in this part some algebraic known results. 

• In Section[31 we explicitly compute E{X, fi, v) with the help of Spitzer's identity. This is Theorem l3.1l 

• Sections m [5] and [S] are devoted to successively inverting the Laplace-Fourier transform with respect 
to /Lt, V and A respectively. More precisely, in Section 01 we perform the inversion with respect 
to jjL] this yields Theorem 14.11 Next, we perform the inversion with respect to v which gives 
Theorems 15.11 and 15.21 Finally, we carry out the inversion with respect to A and the main results of 
this paper are Theorems 16.21 and 16.31 In each section, we examine the particular cases N = 2 (case 
of rescaled Brownian motion), = 3 (case of an asymmetric pseudo-process) and TV = 4 (case of 
the biharmonic pseudo-process). Moreover, our results recover several known formulas concerning 
the marginal distribution of T(t). 

• The final appendix (Section [7]) contains a discussion on Spitzer's identity as well as some technical 
computations. 

2 Settings 

2.1 A first list of settings 

In this part, wc introduce for each integer n a step- process A„ coinciding with the pseudo-process X on 
the times fc/2", fc G N. Fix n e N. Set, for any fc € N, Xk,n = X{k/T') and for any t e [fc/2", (fc + l)/2"), 
X{t) = Xk^n- We can write globally 



Now, we recall from [13] the definitions of tame functions, functions of discrete observations, and admis- 
sible functions associated with the pseudo-process X. They were introduced by Nishioka [T5^ in the case 
= 4. 

Definition 2.1. Fix n G N. A tame function for X is a function of a finite number k of obser- 
vations of the pseudo-process X at times j/2", 1 < J < fc, that is a quantity of the form Fk^n — 
-F'(A(l/2"), . . . , X(fc/2")) for a certain k and a certain bounded B or el function F : K'^ — > C. The 
"expectation" of Fk^n is defined as 



oiX. 




OO 





F{xi, . . . ,Xk)p{l/'2";x - xi) . . .p{l/2'^;xk-i - Xk)dxi . ..dxk- 
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Definition 2.2. Fix n G N. A function of the discrete observations of X at times k/2'^, k > 1, is a 
convergent series of tame functions: Fx^ = X^feLi where Fk^n is a tame function for all k > 1. 
Assuming the series XifcLi |lE(Ffc,n)| convergent, the "expectation" of Fx„ is defined as 

oo 
k=l 

Definition 2.3. An admissible function is a functional Fx of the pseudo-process X which is the limit 
of a sequence {Fx„)neN of functions of discrete observations of X : Fx = lim,i->.oo -F!x„ , such that the 
sequence (E{Fx„))neN is convergent. The "expectation" of Fx is defined as 

E(Fx) = lim E(Fy„)- 
In this paper, we are concerned with the sojourn time of X in [0, +00): 

T{t)= f V+oo)(^(s))ds. 
In order to give a proper meaning to this quantity, we introduce tlie similar object related to X„: 



T„{t)= / l[o.+oo)(^«(s))ds. 
Jo 

For determining the distribution of Tn{t), we compute its 3-parameters Laplace-Fourier transform: 



-At+i/iX„(t)-i/T„(t) 



In Section[3l we prove that the sequence (£'„(A,/i, i'))ner-i is convergent and we compute its limit: 

lim En{X, fi,iy) = E{X, fi,v). 

Formally, E(X,fi,i>) is interpreted as 



-Xt+iiJ,X{t)-iyT{t) 



where the quantity ^-\t+itj.x(t}-uT{t) jg ^^^^ admissible function of X. This computation is performed 
with the aid of Spitzer's identity. This latter concerns the classical random walk. Nevertheless, since it 
hinges on combinatorial arguments, it can be applied to the context of pseudo-processes. We clarify this 
point in Section [3l 

2.2 A second list of settings 

We introduce some algebraic settings. Let 9i, 1 <i < N,he the TV"^ roots of k„ and 

J = {ie {l,...,N} : m^> 0}, K = {i e {1,...,N} : < 0}. 

Of course, the cardinalities of J and K sum to TV: + = N. We state several results related to 
the 0i 's which are proved in [TTJ [T3] . We have the elementary equalities 



N 



jeJ keK i=i 



N 



\keK 



and 



N 



N 



Moreover, from formula (5.10) in |13j. 



keK 



Y^{-lfa,x*^- 

^=0 



(2.1) 



(2.2) 



(2.3) 
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where cre = J2 fci< - </c« ^fci • ■ • Ski- We have by Lemma 11 in [TT] 

ki,...,ktGK 



E^. n 

jeJ ieJ\{j} 



^ jeJ keK 



cos 



2Ar 



2 sin ^ sin f 



if is even, 



if N is odd. 



(2.4) 



Set A, = U^eJ\{J} for j e J, and Bfc = J] 

Vandermonde system: we have 



ieK\{k} e -0k fo'" ^ ^ ^- The Aj's and Bk's solve a 



^A, = 5]B, = 1 

jeJ keK (2.5) 
^ Aj-e'T' = for 1 < m < # J - 1, ^ Bfe6l^ = for 1 < to < #X - 1. 



keK 



Observing that 1/Oj = Oj for j € J, that {Oj^j G J} = {6'j,j G J} and similarly for the ^^'s, k G K, 
formula (2.11) in [13] gives 



UkeA^-Skx) UkeA^-^i^ 



(2.6) 



In particular 



E 



1 

(l-0fcX) 



1 



9j — 9k NBk 



E 

keK 



BkOk 



1 



?fc - 



(2.7) 



Set, for any to S Z, a™ = J2jeJ ^i^T ^^'^ "^"^ ~ "^keK ^^^T ■ have, by formula (2.11) of [O] , 
= {-l)*''-^Y{keK&k. Moreover, = (-l)^^-^!!^^^ ^0 (Efeeif ^fe) • The proof of this 

claim is postponed to Lemma 17.21 in the appendix. We sum up this information and (|2.5p into 



/3m = < 



We also have 



and then 



1 if TO = 0, 

if 1 < TO < #is: - 1, 

(-1)*''-' {I\keK f)^) {TkeK ^fe) if m = + 1, 

K„ if TO = N. 

a-,n = E ^ " ^" E ^J^f "™ = '^ivaAT-m 



(2.8) 



if TO = 0, 



(-1)#^-^ (n,e./ ^.'^ ) ( E,e,7 ) if ™ ^ - 1 





if TO 

if + 1 < TO < A^- 1, 
if TO = A^. 



(2.9) 



In particular, by ()2.ip . 

ao/3o = ol^nPn = 1, a-#KP#K = -1, u-h^kP^k+i = E ^J ' (^i-#kP#k = ^ (2.10) 

With (To = 1, = (DfeeK ^fe) {J2keK ^k) and cr#K = JJ^eK ^k, we also have 

fTo/3o = 1, fT#K-i/3#A- = ^^kP^k+i = (-l)#^^-i ^ ^fe, ^#k/3#a- - (-1)*^-^ (2.11) 
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Concerning the kernel p, we have from Proposition 1 in 

if ^ is even, 

(2.12) 

r(i)cos(^) ^ ' 

if iv is odd. 

Proposition 3 in [TT] states 

poo n J pO 

F{X{t)>0} = j^ P(i;-e)de=^, nX{t)<0} = J p{t;-Od^ = 
and formulas (4.7) and (4.8) in 13i yield, for A > and /i e M, 

poo -At pO ( N/T \ 

foo „ — At ^oo / K/X \ 



N 



(2.13) 



(2.14) 



poo -At i-oo ^ '\ 

Let us introduce, for to < — 1 and x > 0, 
Formula (5.13) in [T3] gives, for < to < — 1 and x > 0, 

i 

e-^^Ij^ra{T]x)AT = e' ^ . (2.16) 



Example 2.1. Case N = 2: we have K2 = +1- This is the case of rescaled Brownian motion. The 
square roots of K2 are — 1, 02 — —1 and then J — {1}, K={2}, Ai — B2 = 1, ao = a_i = 1, /3o = 1, 
= —1. Moreover, 

- / poo poo \ 

The function Ji.o can be simplified. In fact, we have 



Finally, 



h.^'-r:.) = (2,17) 



Example 2.2. Case TV = 3. 

• For K3 = +1, the third roots of K3 are 61 = 1, 62 ^ e*"^', 0^ = e~'^, and the settings read J — {1}, 

= {2, 3}, Ai = 1, B2 = B3 = ^> "0 = a-i = a-2 = 1, /3o = 1, /3-i = -1- Moreover, 

/i,o(r; x) = |- (^^°° e-^e-e^^i _ ^2 ^-.j^-e-'f 

• For K3 = —1, the third roots of K3 are 9i = e't, 6*2 = e^*^, ^3 = —1. The settings read J = {1,2}, 
K = {3}, Ai = A2 - B3 - 1, ao - a_i - 1, /3o = /3-2 = 1, P-i = -1- Moreover, 

/i,i(r;x) = |- (^e-f ^""^e-^^-'^f-^dC-e^t ^e'^^''-< ddj, 
2^ 
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Actually, the three functions Ii^, Ii^i and I2.1 can be expressed by mean of the Airy function Hi defined 

as Hi(z) = i 6"^"+^^ d^ (see, e.g., [H Chap. 10.4]). Indeed, we easily have by a change of variables, 
differentiation and integration by parts, for t > and z € C, 



^e-^^dC- 



(3r)4/3 

" HiY " 



Therefore, 



-fi,o(T;a;) = 
/i,i(t;x) = 
^2,i(t;x) = 



e'eHi 



(3r)2/3 



/3r 



Hi 



(3r)4/3 3r 



e ax 



e-*6Hi - 



e 3 X 



/3r 



2r2/3 
2t2/3 



/3t 



e'eHi' 



/3t 



e-*6Hi' 



^3t 

- 271 



(2.18) 
(2.19) 
(2.20) 



Example 2.3. Case = 4: we have K4 — —1. This is the case of the biharmonic pseudo-process. The 
fourth roots of K4 are — e~'^T , 62 — e*^, 63 — , 64 — e^*T and the notations read in this case 

J = {1,2}, K = {3,4}, Ai^ B3 = A2 = Bi = ^, ao = a_2 = 1, «-! = \/2, /3o = /3-2 = 1, 

= -\/2. Moreover, 



2 g-rr+^c 



/2,i(T;x) = ^^e't re 



2 ^-Tj-'-ia;? 







d^ . 



(2.21) 



3 Evaluation of E{X, fi,^) 

The goal of this section is to evaluate the limit E{X, fi, v) — lim„_yoo En{\, fi, v). We write En{X, fj., v) 
E[F„(A,/x,:y)] with 



Let us rewrite the sojourn time r„(i) as follows: 



-Xt+ifj.Xn(t)~vT„{t) 



n([2"t] + l)/2" 



[2^1 rU + l)/2" H[2"t] + l)/2" 

Tn{t)^22 / l[0,+oo)(^n(s))ds - / l[o^+oo){X„{s)) ds 

j=0 "'i/2" Jt 

= 51 / l[o,+oo)(^j,«)ds- / 

= ^Y1 l[0,+oo)(^J,n) + U — ) l[0,+oo)(^[2"t],rO- 



l[0,+oo)(-'^[2"t],ra) ds 



j=0 



Set To^n = and, for fc > 1, 



1 

Tk,n = ^ X/ l[0,+oo)(-'^i,n)- 



For > and t e [fc/2", (fc + l)/2"), we see that 

Tnit)^Tk,n+ (t- 



k+1 
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With this decomposition at hand, we can begin to compute Fn{\, ^,,v): 

g-At+iAi^„(t)-i/T„(t) 
oo /.(fe+l)/2" 







^-At+v^i,.„-iyTfc,„-^+iy(i|fci-i)l[o. + oc)(Xfc.,i) 



oo „(fe+l)/2" 



E 



fe/2" 



The value of the above integral is 

n(fc+l)/2" 



-At+i/(ifci-t)l[o, + oo)(-f|c,„) = g-A(fc+l)/2' 



fc/2" 



e[A+i^l[n. + co)(Xfc,„)]/2" _ 
^ + '^l[0,+oo)(-'^fc,n) 



Therefore, 



1 _ „-(A+i/)/2" oo 



A + 1/ 

-iy/2 



k=0 

1 - e-V2" 



A 



fc=0 



Before applying the expectation to this last expression, we have to check that it defines a function of 
discrete observations of the pseudo-process X which satisfies the conditions of Definition 12.21 This fact 
is stated in the proposition below. 

Proposition 3.1. Suppose N even and fix an integer n. For any complex A such that 3fJ(A) > and 
any v > 0, the series Er=o e"^''^^^"lE[e'^-^'='"-'''^'='"l[o,+oo)(^fe,«)] and Y^'^^^e-^^'^^¥.[e^f'^'''--''^'^-- 
l(-oo.o)(-'^fe,n)] are absolutely convergent and their sums are given by 



fc^O 

CSO 



fc=0 



e-/2"-5+(A,M,^) 
e-/2"[5-(A,A^,i.)-l] 

g^/2" _ I 



where 



I oo 

5+(A,/x,^.)=exp -^(l 



-vkl2'' 



k=l 



-Afe/2" 



Sn i\ M> i^) = exp ^U-e 



^-XkjT^ 



E[e^^^'=."l(_,,,o)(^fe,n)] 



\ fc=l 
Proof 

• Step 1. First, notice that for any fc > 1, we have 



^x[0,+oo) 



R*^-ix[0,+oo) 



< 



dyi . . . dyfe = ^ 



dxi . . . dx/c 



fe 

n 
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Hence, we derive the following inequality: 

oo oo 



\k/2" 



k=l 



k=l 



1 



1 — pe' 



-SR(A)/2" 



We can easily see that this bound holds true also when the factor l[o,+oo)(-''^fc,n) is replaced by l(-oo,o)(-'^fc,n)- 
This shows that the two series of Proposition 13.11 are finite for A G C such that pe~^^^^/'^"' < 1, that is 
5R(A) > 2" log p. 

• Step 2. For A e C such that 5R(A) > 2"logp, the Spitzer's identity ([7^ (see Lemma O in the 
appendix) gives for the first series of Proposition 13. II 



fc=0 



exp 



SO 

k=l 



.-vkjT 



E[e*^^'=-l[0,+oo)(Xfe,n)] 



(3.1) 



The right-hand side of (I3.ip is an analytic continuation of the Dirichlet series lying in the left-hand side 
of p.ip . which is defined on the half-plane {A € C : 3fJ(A) > 0}. Moreover, for any e > 0, this continuation 
is bounded over the half-plane {A e C : 3fJ(A) > e}. Indeed, we have 



|E[e^''^^."l[o,+^)(Xfe,„)] 



+ 00 



2" 



< 



+ 00 



and then 



exp 



fc=i 



< exp ■ 

V fe=i 



-Afc/2" 
-3?(A)fc/2' 



{e'^''''-'^lO.+oo){Xk,n)] 

-J =exp(-plog(l-e-^(^)/2" 



Therefore, if g?(A) > e. 



exp 



k=l 



,-i/fc/2' 



,-Afc/2" 



E[e*^^''=^"]l[0,+oo)(^fc,n)] 



(1 - e-5R(A)/2")p- 



< 



(l_e-e/2")p- 



This proves that the left-hand side of this last inequality is bounded for 9fi(A) > e. By a lemma of Bohr 
([!]), we deduce that the abscissas of convergence, absolute convergence and boundedness of the Dirichlet 
series X]fe°=o ^~'*''^^^"'^[^*^'^'''"~''"^'° "-'^[o,+oo)(-^fc,n)] are identical. So, this series converges absolutely on 
the half-plane {A G C : 3fi(A) > 0} and (I3.ip holds on this half-plane. A similar conclusion holds for the 
second series of Proposition 13. II The proof is finished. ■ 

Thanks to Proposition 13.1) we see that the functional i^„(A,/i, i^) is a function of the discrete obser- 
vations of X and, by Definition 12.21 its expectation can be computed as follows: 



^_g-(A+.)/2" 6-^/2" -5+ (A, /x,i.) l-e-V2" S-{X,p,iy)-l 



1 



A 



1 



,-A/2" 



e'^/'"(l-e-(^+'-)/^") 
(A + i/)(e''/2" - 1) " A(e''/2" - 1), 
1 _ e-^/^" _ 1 - e-(^+'')/2" 



A(e'^/2" - 1) 



(A + i^)(e'^/2" - 1) 



= !'/2" 



'5,|(A,p,z^). 



(3.2) 



Now, we have to evaluate the limit E{X,fj,,i>) of _E„(A,/i, i^) as n goes toward infinity. It is easy to see 
that this limit exists; see the proof of Theorem 13. II below. Formally, we write E{X,^,,v) — E[_F(A, /i, i^)] 
with 

/•oo 

F{X,p,v) = / e-^t+^^.x(t)-.T(t) 
Jo 

Then, we can say that the functional F{X, p, v) is an admissible function of X in the sense of Definition 
The value of its expectation E{X, p, v) is given in the following theorem. 
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Theorem 3.1. The 3-parameters Laplace- Fourier transform of the couple {T(t),X{t)) is given by 

E(X,ii,iy)^ ; — . (3.3) 

Proof 

It is plain that the term lying within the biggest parentheses in the last equality of p.2p tends to zero as 
n goes towards infinity and that the coefficients lying before S^{\,^i,v) and S^{X,^,i') tend to l/i/. As 
a byproduct, we derive at the limit when n — )■ cxd, 

E{X,^i,:.) = - [S-{X,^l,,y)-S+{X,^i,,.)] (3.4) 

where we set 

5+(A,A.,;.) = 5+(A,Ai,;.) ==exp(^-^"E[e*^^(*)l[o,+o,^ 

S-{X, ^l, v) = Jim ^-(A, /i, v) = exp(^ ^ E[e^^^(*)]l(_oo,o)(^(i))](l - e""*) ^ dt^ • 



We have 



E[e'^^«l[o,+oo)(^(t))](l - e-^')^dt 







E[(e*^^(*)-l)%+^)(X(t))] -—dt- E[(e^''^W-l)l[o,+oo)(^(i))] 



oo -At _ -{\+u)t 

¥{X{t) > 0} dt 



dt / - 1) Pit; -0 dC - / dt / (e^^^ - l) p(t; d^ 



^ "'0 

/•oo -At _ -(A+iy)t 

+ P{X(1)>0}/ dt. 

Jo i 



In view of (|2.13p and p.l4p and using the elementary equality ^ ^t'^ ^ ^ '^^ ~ -^'^S (^a^)' have 



-At 

E[e'^^W]I[o,+oo)(^(i))](l - e-^*) ^dt 



We then deduce the value of S+{X, p., v). By ([^ . 

S'+fA u zy) = TT ~ '^•'^^ = n£ti(yA-iiAJ6i£) 

^ A - 

Similarly, the value of ^^(A, /i, v) is given by 



(3.5) 



feeK 'sf^-ilJ'Sk Y{j(.j{'VxT^-ilJ'Oj)]\^^j^{VX~i^iek) 
Finally, putting p.5p and p.6p into p.4p immediately leads to p.3p . ■ 

Remark 3.1. In the particular case ^ = 0, we get the very simple result: 

/.oo 1 

E{XM= e-^*E[e-''^(*)]dt=-^^- 

Jo A « (A + z^) N 

This is formula (20) of 11 . On the other hand, we can rewrite p.3p as 
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E{x, /X, = ^ ^ — -#T n 



n 



(3.7) 



Actually, this form is more suitable for the inversion of the Laplace-Fourier transform. 

In the three next sections, we progressively invert the 3-parameters Laplace-Fourier transform E{X, fj,, v). 

4 Inverting with respect to ^ 

In this part, we invert E{X^^,v) given by p.7l) with respect to /i. 
Theorem 4.1. We have, for X,v > {), 



e-^* [E(e-''^(*\ X{t) e dx)/dx] At 





kOk 



ekVx-0jVxTl^ 



(4.1) 



Proof 



By (j2.6l) applied to a; = iii/ \/X + v and .t = ipi/ \/A, we have 



n 



i6 



n 



Vx 



n ^ n ' 



Ml 

i 

Va(aT^5] 



AjBkOjOk 



jeJ y"3 



V^+^-it^WkVx-iti) 



Let us write that 



,VAT1^ - iM)(6'fcVA - i^i) 6ifcVA - Sj^TTm V^'jVaT^ - i^l OkVX - ijj. 



1 



Ok^-e^'i/xTTi \Jo 

Therefore, we can rewrite E{X, as 

1 



A N (A -|- I^j N 



e*"" (e~''=^"l(-oo,o](^) +e-'^^^l[o,oo)(a:)) Ax 



kGK 



which is nothing but the Fourier transform with respect to of the right-hand side of (|4.ip . 



Remark 4.1. • By integrating (|4.ip on (— oo,0], we obtain 



E(e-''^(*\ X(t) < 0)dt 



A'^(A + 

Using ()2.6p applied to x = Ok'VX/ \fX+ly and (|2.3p . we see that 



kVx~6,VxT^ 
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je.J — ffc 



n/ a 



yiPk 



iei< 



A 



X + v 



n 

#K 



X + v 



This entails that 



lerx^{x + .)^Y.^-iY--^~^r 
^^x^{x + .)^Y.^-iY-^ei 



-I 



X + v 



X + u\ " 



e-^*E(e-''^(*), X{t) <0)dt=-J2 BkJ^i-^Y^i^^i 



keK e=o 



X + iy 



1 „ / X + v 



X 



By (|2.8p . we know that ah the Pi, 1 < (- < #K — 1 vanish and it remains, with (I2.1ip . 



-^*E(e-''^(*), Xit) <0)dt - 



X + v 



X + v 



X 



- 1 



We retrieve (30) of [H]. 
• Likewise, we have 



(4.2) 



(4.3) 



;-^*E(e-''^W,X(i)<0)dt=- 



X + v 



which coincides with (29) of 

• Adding formulas (|4.3I) and (14.41) we obtain 



e-^*E(e-^^(*))dt = 



1 



A + J/ 
A 



A 



X + v 



(A + v)—Ti — - A 
vX N [A + V) N 



(4.4) 



An [X + V) n 



This is formula (10) of [TT which has already been pointed out in Remark |3. II Another way of checking 
this formula consists of integrating (14. ip with respect to x directly on R. Indeed, 



1 



, #.7-1 

A « (A + z^) N 



= E 



AjBkOk 



-—y 



AjBkdj 



keK 



hVx - e,^Vx + D Vx fl OkVx - em 



3&J 
keK 



-J y 



x^{x + v)^ VaTT; (0fe^ - 0j VXT^) Va (0fc^ - e.VATT;) ^ 



1 



AjBk 



x'^iX + v)'^ f^j V^+'^t 
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By ([231), we have j&J ^j^k = Y^^eJ T^keK -^fe = 1 and then 

keK 



e"^*E(e-''^(*))dt 



An (A + i^j n 



Remark 4.2. By replacing a: by into (|4.ip and by using (|4.2p . we get 



fee A' f=o ^ 
= 2^(-l)V,/3,+i ( ^ 

In this sum, all the /3£+i, < ^ < — 2, vanish and it remains, with (|2.1ip 



e 1^ e 



X{t) e dx)/dx 



dt 

x=0 



keK \ / Vjej / 

We retrieve formula (26) of [Tl] . 

Remark 4.3. For z^ = 0, formula (|4T1) yields with (|2?7l) the A-potential of X: 



e-''[P{X{t)edx}/dx]dt^} ^ ^^.^ 



Kg-". VAX if2.>0, 



1 ^-Sj v^x 

We retrieve (12) of pT]. 

Example 4.1. For N — 2, formula ()4.ip gives, with the numerical values of Example 12.11 



/>oo 

/ e-^* [E(e-^'^(*), X(t) e da;) /da;] dt 
Jo 



|\/A + v 



VA+1^ 



/ A+L' a; 



3^=" ifx<0. 



This is formula 1.4.5, p. 129, of 0. 



Example 4.2. For = 3, we have two cases to consider. Formula ()4.ip yields, with the numerical values 
of Example 12. 2[ in the case K3 = 1, 
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e-^* [E(e-''^(*\ X{t) e dx)/dx] dt 



•/X+ux 



= < 



A2/3 + ^A(A + u) + {\ + J/)2/3 

e^- V3 ^cos(^a;)~(2^AT^+^)sin(^a;) 
A2/3 + ^A(A + i.) + (A + vf/^ 



if X > 0, 



if X < 0, 



and in the case K3 = —1, 

[£(6-"^^*^ X{t) £ dx)/dx] dt 
' e-^- V3^AT^cos(^^S^a;) + (^AT^+2^)sin(^^|S^a;) 



A2/3 + ^A(A + j.) + (A + !y)2/3 



[ A2/3 + ^A(A + I/) + (A + J/)2/3 

Example 4.3. For iV = 4, formula (|4.ip supplies, with the numerical values of Example 
e-^* [E(e-"^(*), X{t) e dx)/dx] di 



if X > 0, 
if a; < 0. 



^aT^(x/a + VaT7^)(^+ ^aTT^) L 



V A + cos ^ — a; + V A sm — — x 



r- 

V2e ■ 



I ^(x/A + x/AT1^)(^+ ^ATT^) 



"i^A cos f — v^A + sin ( x 

V2 J VV2 



if X > 0, 



if a; < 0. 



Remark 4.4. Using quite analogous computations to those of Remark 14. 2[ we could derive another 
expression for formula (14. ip . Actually it will not be used for the inversion. 



-At 



[E(e-''^(*\ X{t) e dx)/dx] dt 



E^A n 



if X > 0, 



je.J \te.i\{j} ' ^ 



n^A n 

k£K \ieK\{k} 



(4.5) 



A 



7i — Ok 



Formula ()4.5p looks like formula (24) in [TT]. Nevertheless, (|4.5p involves the distribution of {T{t),X{t)) 
when the pseudo-process starts at zero while (24) of [TT] involves the density of {T{t),X{t)) evaluated at 
the extremity X{t) = when the starting point is x. Actually, both formulas are identical by invoking 
the duality upon changing x into —x, but they were obtained through different approaches. 

5 Inverting with respect to u 

In this section, we carry out the inversion with respect to the parameter u. The cases x < and x > 
lead to results which are not quite analogous. This is due to the asymmetry of our problem. So, we split 
our analysis into two subsections related to the cases x < and x > 0. 



5.1 The case x < 

Theorem 5.1. The Laplace transform with respect to t of the density of the couple {T{t),X{t)) is given, 
when X < 0, by 
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[P{T{t) e ds,X{t) e dx}/{dsdx)]dt 



^ m— 



(5.1) 



Proof 

Recall (|4.ip in the case a; < 0: 



/>C30 

^0 

= T*Eir;^—T*ni Yl ^^^^ ( ■ 



We have to invert with respect to v the quantity 



By using the following elementary equality, which is valid for a > 0, 
1 1 



(A + i.)" r(a) io 



r(a) 



we obtain, for |/3| < \/A + 
1 1 



ds, 



r=0 V Af ^ 



The sum lying in the last displayed equality can be expressed by means of the Mittag-LefHer function 
(see Chap, xviii]): Ea,b{0 = EZo rS+b)- ^hen, 

1 



Next, we write 



A, 



J e~""' i^sir-^e-^'Ei.^^{P^)jds. 



where 



i6J 



jGJ r=0 



J U_. ] ^A^\oJ r(^) ^0 ^ ; r(i±i) 



ds, 



A, \ (A,s): 



(5.2) 



(5.3) 



AT 



When performing the euclidian division of r by N, we can write r as r = £N + m with i > and 
< m < iV - 1. With this, we have 9 J'' = {Of )-'^ Oj"' = OJ"" and 01 ^k^^O^. Then, 

^I- E # - E ^ = 



Hence, since by p.9p the a_m, + 1 < m < A^, vanish, 



E^.^^,^ =EE 



(As) 



£=0 m=0 
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and (j5.3|) becomes 

A r°° ( *^ \ 

As a result, by introducing a convolution product, we obtain 



e-^* [E(e-"^(*), e dx)/dx] di 



X 



/ '"''M / X ^'"'^'"^ E «->»^rA^(g - a)'^-'E, ,^{\{s - a)) da ) As. 

■^0 ^ V W j m=0 / 

By removing the Laplace transforms with respect to the parameter v of each member of the foregoing 
equality, we extract 

/>C30 

/ e^^* [F{T(t) e ds, e da;}/(dsda;)] At 
Jo 

^ E a_™At [ ^ B.^^+ie-^^-v^-) - a)"^-^^!,^ (A(s - a)) da. 



A « 



The integral lying on the right-hand side of the previous equality can be evaluated as follows: 



^(s-a) « i£;^^^(A(s-a))da= / — ^(s-a) « i;^^^_^^^da 



-1 



— . „ = S " n.+#j (As) 

Y p I m+#J \ 1. w ^ ^ 

^=0 AT 



from which we deduce (j5.ip . 



Remark 5.1. Let us integrate formula (|5.ip with respect to a; on (— oo,0]. This gives 

poo -^Xs #K 

/ e-^* [P{T(t) e ds, < 0}/ds]di = ^ E a-m/3™(As)^£;i ™+#i(As). 

(-^s) " m=0 

In view of (EH) and (PJU)) . since i;i,i(As) = e^", 

/ e-^* [P{T(i) e ds, X{t) < 0}/ds]dt = (aoPoE^ #i (As) + a_#if/3#K(As) — £i,i(As) 

Jo (As) N ^ ■ ™ 



(As)^ 

l^E,^^{Xs)~l. (5.4) 



(As) 

We can rewrite E-^ #j(As) as an integral by using Lemma [731 We obtain that 



N 



/ e"^* [P{T{t) e ds, X{t) < 0}/ds] dt 
Jo 



1 La. , r*'-^„-A.*dt 



e 



(As)^ (As)^-i V Jo i+1 



dt^ ''2k ' e' / ^ e-^' dt 

Jo ^ + 1 TTS N Jo ^ ~^ ^ 



TTS « 



dt. 
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From this, we extract, for < s < t, 



sin(^^7r) /t—s 



p{rm e ds, < o}/ds = — ' ^ M . (5.5) 

T:t \ s J 

We retrieve Theorem 14 of [TT]. By integrating (|5.5p with respect to s, we get 



which simphfies to V{X{t) < 0} = #A77V. We retrieve (11) of [H]. 
Remark 5.2. An ahernative expression for formula (|5.ip is for x < 



oo 

At 



g^At [p{y(t) £ ds,X(t) e da;}/(dsda;)]dt 

1^ ^ _^ f^^) e-«^"^^ (5.6) 



In effect, by (fSTTj) . 

/ e-^* P{T(t) e ds,X{t) e da;}/(dsda;)]dt 
Jo 



-As #^ ^'\„^^+■^ 



m+#,J- 



g Wfc VAa: 



1=0 m=o \ J / i i^t -I- ^ 



In the last displayed equality, we have extended the sum with respect to m to the range < m < TV — 1 
because, by (|2.9p . the a_m, + 1 < m < N — 1, vanish. Let us introduce the index r = £N + m. Since 

j;) ,wehave 



/ e-'' mm E ds,X{t) e d:r}/(dsd.T)]d< = J] ^.^fe^fe E ^^77T#7r 

which coincide with ([5~ 

Example 5.1. Case N — 2. Suppose a: < 0. The first expression (|5.ip reads 

.,-As 



/ e"^* [P{T(i) e ds,X(i) e da;}/(dsda;)]dt = i(As) - VAs (As) ) 



while the second expression (15.61) reads 

foo „ — As 



/•oo p-As 

/ e"^* [P{r(i) e ds,X(t) e da;}/(dsda;)] dt = i (- Vx^) e^^^. 



From Lemma 17.51 '^6 have 

£;i.i(-^/A;) = £;i,i(As)-VAie^^ (5.7) 
which proves the coincidence of both formulas. Moreover, from Lemma |7.5[ for a; < 0, 

/ e-^* [P{T(t) € ds,X(t) e da;}/(dsda;)]dt = ( - \/A Erfc(\/Al)^ e^^. 

io Vv^ / 

We retrieve formula 1.4.6, page 129, of [B]. 
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Example 5.2. Case N = 3. We have for a; < 0, when K3 = — 1: 



-As 



2(As) - 



and when K3 = 1: 



/ e"^* e ds, e da;}/(dsda;)] di 

"'0 



73 



V3 cos(' ^^'^ V^£;i,|(Ag) - (As)2/3< 



As 



sm(^^^^^^(^£;i^. (As) + (As)2/3e^^ - 2E,^^{X.s)) 



Example 5.3. Case = 4. We have, for x >0, 



,-At 



"{Tit) e ds, X{t) e dx}/(dsda;)]dt 



cos( ^S^V 3 (As) - VAie^^) +sinf^^V^£;i 3 (As) - 1 (As) 



V V2 



V n/2 



5.2 The case x>0 

Theorem 5.2. T/ie Laplace transform with respect to t of the density of the couple {T{t),X{t)) is given, 
when a; > 0, by 



poo 

/ e'^^[¥{T{t) £ ds, X{t) e da;}/(dsda;)]dt 
"'0 



-As 



^Y^A^BkOk / a^-'E^J^yX^]l,^#j^,{s-a-x)da 



where the function Ij^^j^i is defined by i2.15\) . 
Proof 

Recall (|4.ip in the case a; > 0: 

/>oo 

/ e-^* [E{e-'^^'\ X{t) e Ax) /Ax] At 
Jo 



A « (A + i^) « ^.gj Vfcg^ 



BkOk 



(5.8) 



We have to invert the quantity 



(A-H^)^^(«/A+i;_|i.^) 
^ ^ (sTr-^e-^'Ei. !. (/3 Vi) ) ds 

(e-^%-,#,7_i(s;x))ds 



with respect to v. Recalling ()5.2p and ()2.16p . 



vatt;-/? jo 

g-e^VA+I^a; 
{\ + V) N Jo 



we get by convolution 
(AT^^)*^^^^XT^"^]^Wy 



E.. VA^ X e-^^'^--)/,- #^_i(s - a; a;) da ds 
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This immediately yields ([S 



Remark 5.3. Noticing that 



\ -J ^ r=0 J ^ V W J f=o m=0 i I* + W j „i=o J 

and reminding that, from (I2.8p . the 1 < m < ^j^A' — 1, vanish, we can rewrite (I5.8P in the following 
form. For x > 0, 



CO 



e 



[P{T(t) e ds, € dx}/(dsdx)] dt 







E E^^^rMA^^ ra-^-x-^(A-) E#^*-i('^--;-) 



= — e 

N - - 

= -e-^' E ^rnX"^ I a^-^Ei,^{X<j)^,n{s~-<j;x)da (5.9) 

m=H=K 







with $„(r; a;) = ^^.^^ -^^ /j- #j_i(t; x). 

J 

Remark 5.4. Let us integrate (|5.9p with respect to x on [0,oo). We first compute 



r(^) sin(^.) 



TTT N 1 r JV 



Then, with the aid of (H^H) and (PtU)) . we get 



e 



[P{r(i) e ds, > 0}/ds] dt 







a^^K I3#K -^E^ #K{\a)da + a^N I3n X « / ■ —da 

Jo (s ^ a) N ' " Jo (s ~ a) N > 



r(^) V Jo (s-a)^ ' Jo {s~ay 

■da- X« ^ da 



V'1v"j V"^" a N i^s — a) N Jo [s — a) N j 



= e > -^r TiTTT / do- = 1 -ry- / a— " da. 

By Lemma l7.3[ this simplifies into 

g-At [f>{T{t) e ds, X{t) > 0}/ds] dt = 1 - {Xs)^e'^'E^ ,tf+i{Xs). (5.10) 
Now, using Lemma 17.41 we derive another representation for the foregoing Laplace transform; 



sm 



g-At [p{T(i) g ds, X(t) > 0}/ds] di 1 - e--'^ e^'^ / e"-'^* dt 

\ TT 7o t + 1 ^ 

sin(^7r) #,/ e-^(*+^) , 

V w ^ „— / dt 



(t + S)tN 



sin(^7r) ^ 



s « / — dt. 

s t{t ~ s)~ 
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As a result, we derive 



"{Tit) e ds,X{t) > 0}/ds = 



TTt 



t - S 



0<s<t. 



(5.11) 



This is formula (11) of [TT]. 

Remark 5.5. If we add and ([ETU)) . we find 



"{Tit) e ds}/ds] dt = 



E 



E 



-As 



^r(^ + %^) ^r(^ + %^ + i); r(^)(A.)^ 



It is easy to invert this Laplace transform. Indeed, 



-A.S 



-As ^oo 



X iv r(^) Jo 



t^-^e-^' dt = 



#K ■ 
N , 



(t- s) " 'e 



^^-^'dt 



n \ N 



sm 



7 



oo g-At 



(t-sy 



■ dt. 



This implies that 



"{Tit) e ds}/ds = 



sm 



s N (r — sj N 



which can be also obtained by adding directly (|5.5p and ()5.1ip . Thus, we retrieve the famous counterpart 
to the Paul Levy's arc-sine law stated in [11 (Corollary 9). 

Remark 5.6. For .t = 0, using formula ([STTjl which is valid for x < 0, we get, by (|2^ . (|2^ and ()2.10p . 



'{T{t) e ds, X{t) e d.x}/(dsda;) 



dt 



-As 



, *K-1 #K 
A N S JV 



Ea_™ /3™+i (As) " E, r,.+#j (As) 



-As 



#/c-i #^ (As) " " £:i i_^(As) +a_#;f/3#K+i(As)"" i;i4(As) 



A « s jv 

„-As 



A N s « 



5]0,(As)*^i?i,i_^(As) + eu{\s)^e^^ 



E^. 



-As 



^i.i-i(As) 



(5.12) 



On the other hand, with formula (|5.9p which is vahd for x > 0, 

/•oo 

/ e-^*P{r(t) e ds, X(t) e da;}/(dsda;) dt 
Jo ^=0 

PmX"^ / fT^-i^i,^(Aa)$,„(s-a;0)d(T 

„-jii.' "'0 



(5.13) 



m=#K 



with 



\jeJ 1 / 



r(^)Bin(^.) 
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In view of and ([2l0l) . we have 

/ e^^* P{T{t) e ds, X{t) e dx}/(ds dx) dt 



-As 



(s — a) 



■ #ii (Act) da 



(S — (Tj 



■ -£'1 (Act) dcr 



OO 

f=0 



-As 



,-As 



E 



-° AT ' Af 



.0 + 



N 



,-As 



-'/TTgAs 



Thus, we have checked that the two different formulas (|5.12l) and (I5.13P lead to the same result. 

Example 5.4. Case N = 2. Suppose x > 0. Formula (|5.8p reads, with the numerical values of Exam- 
ple 

poo pS 1 

/ e-^'[V{T{t) €ds,X{t) €dx}/{dsdx)]dt = e~^'' -= Ei i{- VX^) h^Q^s - (t;x) da 
Jo Jo vo" ^'^ 

while formula ()5.9p gives, because of$i=$2=/i,o and ([H 



-At 



^{T(t) e ds, X(t) e da;}/(dsda;)] dt 



-As 



.-As 



/ 1 (Act) $i(s - ct; a;) do- - \/A / i(Act) $2(5 - cr; x) dcr ) 

Jo Jo ' J 

I -^(Ej^i (Act) - Vx^e^") h o{s - cr; x) dcr 
Jo Vcr ^ / ' 

/ — = i?i 1 f— VXa) h ois — cr; x) da. 
Jo Vo^ 



We have checked that the two different representations (j5.8p and (|5.9p lead to the same result. Let us 
pursue the computations. In view of p.l7p and Lemma [7.5) we get 



g-At [p{2^(t) g ds, € da;}/(dsda;)] di 



-As /-s 



g 4(s-<t) 



2%/^ Jo V^(s-a)3/2 



£■1,1 (- \/Act) dcr 



-As ps 



^1 El 1 (- v/A(s ~a))da 

TT Jo a^l^^s - a ^ ^ 



— As i-s 



e 4" 



1 



2%/7i^ 7o a^/'^^/s - a VV^ 



- VA(s-ct) e^("-") Erfc(VA(s-a)) ) da 



-As 



2V^ 



1 



e 



dcr - 



2%/A r e^(^-'^)- 



TT Jo a^/'^^/s — a 0r 7q 



73/2 



.e-«'dn dcr 



The first integral in the last displayed equality writes, with the change of variable a = /t 

2 

1 



:dcr 



(t3/2VI^ s3/2 ^ 



£43'! g 4 

= dr 



3/2 



e , 2^/7? _^i2_ 
-j=— dr = e 4= ^ 
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and then 



/ e"^* [V{T{t) e ds,X{t) e da;}/(dsda;)]dt 



1 a;%/A 

e 4= 



The computation of the integral lying on the right-hand side of the foregoing equality being cumbersome 
is postponed to Lemma 17.71 in the appendix. The final result is, for a; > 0, 

/ e"^* [P{T(t) e ds.Xlt) e da;}/(dsda;)]dt = VAe^^^Erfc — ^ + VaI . 

Jo VTTS \2Vs J 

This is formula 1.4.6, page 129, of [6]. 

Example 5.5. Case = 3. For x > 0, (|5.8p supplies with the numerical values of Example 12.21 when 
K3 = -1, 

/ e"^* [P{T(t) G ds, € da;}/(dsda;)]dt 

= e- / cr-2/3i;i 1 ^)/ll(s-cr; x) da 

v3 V Jo 

+ e-'s' ^ a-^^^Ei i (- /2,i(s - ct; x) dcr 



and when K3 = 1, 



00 



Jo 



g-At [p{2^(t) (= ds, e da;}/(dsda;)] dt 

- ^ cr"^/^£;i , 1 (e*^ ^) /i,o(s - cr; x) da^ . 

The functions Ii^o, and /2,i above are respectively given by (|2.18p . (|2.19p and (|2.20p . 
Example 5.6. Case iV = 4. For a; > 0, (|5.8p supplies, with the numerical values of Example 



/ e"^* [P{T(i) e ds, X{t) e dx}/(dsda;)] dt 
Jo 



2^ 



+ e^'T ^ a-^/*Ei i (- ^) /2,i (s - cr; x) da^ . 
The functions Ii i and /2,i above are given by (|2.2ip . 

6 Inverting with respect to A 

In this section, we perform the last inversion in F{X, fi, v) in order to derive the distribution of the couple 
{T{t),X{t)). As in the previous section, we treat separately the two cases x < and x > 0. 
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6.1 The case x <0 

Theorem 6.1. The distribution of the couple (T{t),X(t)) is given, for x < 0, by 

P{T{t) e ds,X{t) e dx}/dsdx 

AT ' 

= -^E"-™-"^/ r+^"e-(-^)«"^™(xOi?i^(-.?^)cie (6.1) 

where 
Proof 

Assume x <0. Recalling (|5.1I) . we have 



g-At [f>{T{t) e ds, X{t) e da;}/(dsda;)]dt 

-As 

m=0 f=0 ^" N ) keK 

= - E E "-'"F^TTSZY E A^""^e---^"v^^ (6.2) 

^=0 m=0 -L (^t + ^ ) 

We need to invert the quantity A^^ * ^ g-As-Sfc \/Aa; for £ > Q and < to < with respect to A. 
We intend to use (|2.16p which is valid for < to < iV — 1. Actually (|2.16p holds true also for to < 0; 
the proof of this claim is postponed to Lemma 17.81 in the appendix. As a byproduct, for any £ > and 

< TO < #X, 

^,+ n,^^.,s-e.'VX.^^-Xs / e-^"4,#K-.A.-™-i(^;2:)d^ 



e Ik,i^K-(N-m-i{t - s;x)dt. (6.3) 

Then, by putting (|6.3p into (|6.2p and next by eliminating the Laplace transform with respect to A, we 
extract 

P{r(s) e ds, X{t) € da:}/(dsd2;) 

= -E E "-"'p/. m+#7 Y E 5fc^'r^^4,#if-£W-m-i(i-s;a;) 

JVl \ - \ ^ s 



IS L X - X - 

X 



/ /"OO ^ 



fee a: 

m=0 fee/f 



; #g-,„-i , 
X e « 



^ J: 
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AT- *^ 

m=0 keK 



X e « 



oo 



The proof of (|6.ip is established. ■ 

Remark 6.1. Let us integrate (|6.ip with respect to x on (— oo,0]. We first compute, by using p.Sp . 



. keK 



if 1 < m < 



We then obtain 



{r(0 G d., < o}/d. = ^'^"^ e-(*-)«" iJ, (-.e^) dc 

ITS « Jo ' " 



TTS « (t — Sj N ^^jj \ t i 

We retrieve (I5.5p . 

Remark 6.2. Let us evaluate P{T(t) e ds,X(t) e da;}/(dsda;) at .t = 0. For < m < #K, 

IC„M = e-'^^^- E B.er' - e'^^^- E = -2»sin( ^^'^^'^ tt) /3„,+i. 

Observing that sin( ^-^]y'"~ V ) = if m = - 1, in view of (gH), (Hill) and (PIU)) . we get 

F{T{t) e ds, X{t) G dx}/ds 

a;=0 



N 






sin 






N 






sin 







Thanks to (|2^ and ([2TT21) . we see that 

P{r(t) e ds, e dx}/ds ^-p{t;0) 

x—O t 

and we deduce 

P{T{t) e ds\Xit) = 0}/ds = liMlM, 
that is, (T{t)\X{t) = 0) has the uniform law on (0,i). This is Theorem 2.13 of [TT]- 
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6.2 The case x>0 

The case a; > can be related to the case x < by using the duaUty. Let us introduce the dual process 
(X(*)(>o of {Xt)t>o defined as = —Xt for any t >0. It is known that (see [II]): 

• If is even, the processes X and X* are identical in distribution (because of the symmetry of the 
heat kernel p): X* = X; 

• If is odd, we have the equalities in distribution (X+)* = X~ and {X~)* = X+ where X^ is the 
pseudo-process associated with = +1 and X^ the one associated with k^^ = —1. 

When N is even, we have {—6j,j E J} = {6k, k E K}. In this case, for any j E J, there exists a unique 
k E K such that 9j = —9k and then 

^^■= n = n r^-T^= n 0-^=^^ 

and 

When N is odd, we distinguish the roots of in the = +1 and = —1: 

• For K„ = +1, let 6I+, 1 < i < TV, denote the roots of 1 and set J+ = {i e {1, . . . , iV} : ^(9+) > 0} 
and K+ ^ {i E {1,...,N} : 5R(6l+) < 0}; 

• For K„ = -1, let 9:[,1 <i<N, denote the roots of -1 and set J- = {i E {I, . . . , N} : ^(9^) > 0} 
and K- = {iE {1,...,N} : ^(9-) < 0}. 

We have {9-,i E J-} = {-9+,k E K+} and {9,:,k E R-} ^ {-O+J E J+}. In this case, for any 
j E J~ , there exists a unique k E such that 9J = —9^ and then 

= n e- -9- " n 1 Q+ = n 0+^5/+ = 

and similarly = S^T- Moreover, we have 

and similarly q:+ = (—1)™/?^. 

Now, concerning the connection between sojourn time and duality, we have the following fact. Set 



T{t)= ]l(o,+oo)(^("))du and T*{t)^ l[o,+oo)(X*(u)) du. 
Ja Jo 

Since Spitzer's identity holds true interchanging the closed interval [0, +00) and the open interval (0, +cxd), 
it is easy to see that T{t) and T{t) have the same distribution. On the other hand, we have 



Tit)= / ]l(o,+oo)(^("))du= / l(_oo.o)(^*(«))d?/= / [l-l[o.+oo){X*{u))]du = t-T*{t). 

Jo Jo JQ 

We then deduce that T{t) and t — T*{t) have the same distribution. Consequently, we can state the 
lemma below. 

Lemma 6.1. The following identity holds: 

¥{T{t) E ds,X{t) E dx}/{dsdx) ^F{T*{t) E d{t-s),X*{t) E d(-a;)}/(ds dx). 
As a result, the following result ensues. 
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Theorem 6.2. Assume N is even. The distribution of (T{t), X(t)) is given, for x > 0, by 

F{T{t) e ds,X{t) e dx}/(dsda;) 

AT ' r'^ 
= 1^ ^ /3_„ (i - / r^*^ e-«" J^ixO E^^^{-{t - s)e) (6.4) 

^ m=0 ™ 

where 
Proof 

When N is even, we know that X* is identical in distribution to X and (T*(t),X*(i)) is then distributed 
hke {T(t),X{t)). Thus, by (HI]) and LemmaO for a; > 0, 

P{T(t) e ds, e da;}/(dsda;) 

= P{r(t) e d(i-s), X(t) e d(-a;)}/(dsda;) 

= E -^)"^ / r+*"e-^«"/C„(-xOi?i::^(-(i-.s)e^)de 

The discussion preceding Lemma 16.11 shows that 

We see that ICm{z) = (— l)™^^i7m(— ^) where the function J',n is written in Theorem 16.21 Finahy, by 
replacing a-m by (— I)™/?-™ and #J, #_ftr by i^K, #J respectively (which actually coincide since N is 
even), (|6.4p ensues. ■ 

If is odd, although the results are not justified, similar formulas can be stated. We find it interesting 
to produce them here. We set T='=(i) = l^Q^^ao)iX'^ (u)) du. 

Theorem 6.3. Suppose that N is odd. The distribution of {T^ (t) , X^ (t)) is given, for x >0, by 
F{T+{t) e ds, X+{t) e dx}/(dsdx) 



2tt 

m— 



^/3+„(t-s)=^ / r^*^^ e-^" J^ jxQ (-{t s)e) de (6.5) 



where 



J-+(z) = e"' *'^ ^ A+(g+)"'+ie-'^/'^'^^ - e' *'^«"'" " ^ A+(0+)"+ie-''-' ^ 
Proof 

When TV is odd, we know that {X+)* = X- and then [[T+)* {t),{X+)* {t)) = [T~ [t) , X^ {t)) . Thus, 



by (|6.ip and Lemma I6TT1 for x > 0, 

P{T+(t) e ds, X+(t) e dx}/(dsdx) 

= f{T-{t) e d(i-s), X-(t) e d(-a;)}/(dsda;) 

= E "=™(^ - ^)"^ / ^™(--0 (-(i - s)^^) de 



m=0 

where 

/C- (z) = e-*^^^'^ ^ B,7(0,-)™+ie-'^^-^'-^-e'^^^^- ^ i?,-(0,T)'"+ie-^-^'*^ 



As in the proof of Theorem l6.2[ we can write JC'^^{z) = (^ l)™^^i7m (^ 2^) where the function is defined 
in Theorem l6.3l Finally, by replacing a~ by (— 1)™/?+ and #J^, by #1^+, respectively, (|6.5p 

ensues. ■ 
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Formula (|6.5p involves only quantities with associated '+' signs. We have a similar formula for X~ 
by changing all '+' into '— '. So, we can remove these signs in order to get a unified formula (this is (|6.4I) ) 
which is valid for even N and, at least formally, for odd N without sign. 

Remark 6.3. Let us integrate (|6.4p with respect to x on [0, oo). We first calculate, recalling that 
J7m(z) = (— l)'"+^/Cm(— and referring to Remark iG.l) 







if 1 < m < # J, 



Then, 



J-oo — r sin —t^tt if m = 0. 



mt) e d., X{t) > 0}/d. = r C*^-' e-«" E, ^{^{t^ s)^^) 

TT t - s]^ Jo ' " 



TT{t — s) N ^0 

7Vsin(^7r)^(-(t-s)) 



t. /-oo 

Ni+#K-l-sf 



sm(^7rj / t-s\ sm(^7r 



7rt \ t — s 



TTS « — S) « ^^jj 

We retrieve (IS.lip . 
6.3 Examples 

In this part, we write out the distribution of the couple {T{t),X{t)) in the cases N = 2, N = 3 and 
N = A. 

Example 6.1. Case N = 2. Using the numerical values of Example 12. 1[ formula (|6.1I) yields, for x < 0, 
F{T{t) e ds,X{t) £ da;}/(dsda;) = - | ^ 



/>oo 

Jo 

ree-^'-^^i'lC,{xOE,4-se)d^' 
Jo J 



with 

On the other hand, we learn from Lemma [7751 that, for ^ > 0, 



1 



Therefore, 

/•CXD 

P{T(t) e ds,X(t) e d4/(dsdx) = , ' ^ / i (e"«-(*-^)«' - e-"«^(*^^)«' ) d^ 



TT-^'WS Jo 



'e"«-*«%e-"«-*«Mde 



Jo ^ 



: sin(2;0 e-*«' (^J^ ^ e"^ du^ d^ 



^e-«-(*-^)«'df + - / e'e"«-*«'de 
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Let us compute two intermediate integrals: 



2{t-s) 



2(i-s)3/2 



„2 poo 



4*2 



dC 



_1 p 4t 1 

2*3/2^ r 2t 



We deduce the following representation: for x < 0, 

4^ 



"{Tit) e ds,X(t) e da;}/(dsda;) = 



.3/2 



1 



27^^3/2 V 2i 



For a; > 0, IpA^ gives 



1 / 1 /"^^ 

HT(i)eds,X(t)edx}/(dsdx) = - -_ / Ce-^«' JoK)i?i4(-(t-s)^')dC 



TT \y/t — S JQ 
poo 



with 

As previously, 
P{T(t) e ds,X{t) e da;}/(dsdx) 



oo 
oo 



e2e"«-*« d^ 



^ sin(a;0 { du ) d^ 

7r3/2^t - s Jo VJo 



7r3/2 



1- I e 



Actually, the density of {T{t), X{t)) related to (rescaled) Brownian motion is well-known under another 
form. For instance, in [6] pages 129-131, we find that 



' e Ay if X < 0, 



"•{Tit) e ds,X{t) e da;}/(dsda;) = 



/o 47rs3/2(i_ s)3/2 



yiy + x) 



/o 47rs3/2(t_ s)3/2 



e ■!= dy if X > 0. 



The coincidence of our representation and that of ^ can be checked by using Lemma [77^ and Lemma FZ.IOI 
in the appendix. 

Example 6.2. Case iV = 3. 

• Suppose K3 = 1. Using s^"^) = e^^^^ and the values of Example 12.21 (|6.ip writes, for a; < 0, 
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F{T(t) € ds,X{t) e dx}/{dsdx) 
2tt V 



s 



/ ee-(*-^)«'^o(xC)£;i..(-.e')de 
Jo 



-1/3 



/"OO pOO 



where 



^o(^) = -iVSICoiz) = - e--/2 (^cos ^ + x/s sin 

^i(^) = -tVSICiiz) = -e^ + (^cos ^ - VS sin ^ 

^2(^) = -i\/3^2(^) = e' + 2e-"/2 cos 
For X > 0, (I6.5P gives 

P{r(i) e ds,x(t) e dx}/(dsdx) 

where 

^/3z 



:fo(z) = iJo(;2) = 2e-^/2 



sm ■ 



2 ' 

J,iz) = -^J^{z)^ (^VS cos ^ - sin ^) 



• Suppose K3 = —1. Likewise, for a; < 0, 



V{T(t) e ds,X{t) e dx}/(dsdx) 
3/1 
^ 2^ 



1 /"OO /"OO 



where 

V3z 



JCo{z) = -iJCoiz) = -2e^/2 
ti{z) = -ilCiiz) = 



a/S z . \/3 
cos — h sm — - — 



For X > 0, 



"{Tit) € ds,X{t) e dx}/(dsdx) 
2tt 



where 



Jo(^) = ^V3Mz) = e-^ - e^/^ (^cos ^ - VS sin , 
J,(z) ^ ^zVsMz) = -e-^ + e^/2 (^cos ^ + Vs sin ^) 



J2{z)^iV3J2iz)^e-' + 2e''^^cos 
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V3i 



Example 6.3. Case N = 4. Referring to Example 12.31 formula (|6.ip writes, for a; < 0, 
P{T(t) e ds,X{t) G dx}/{dsdx) 

where 

1Cq{z) = —ilCo{z) — — cosz — sinz, 
/Ci (z) = —i ICi (z) = — + cos z — sin z, 
/C2(z) = — i /C2(z) = + cos z + sin z. 

For a: > 0, MM reads 



F{T{t) e ds,X(i) e dx}/(dsdx) 
2 / 1 



/ 1 



V2 



where 



j7o(z) = i JTo{z) ~ e^^ — cos z + sin z. 
J'l (z) = — i jTi (z) = — e^^ + cos z + sin z, 
^2(2) = i ^72 (2) = e^^ + cosz - sinz. 

7 Appendix 

Lemma 7.1 (Spitzer). Let {^k)k>i be a sequence of independent identically distributed random variables 
and set Xq — and Tq = and, for any k > 1, 



Then, for /i G K, > and \z\ < I, 



^E[e^^^'"-''^^]z'= = cxp| 

k=0 



- (.k, ?"fc = l[0.+oo)(^fc) 



00 



,iM^fc-''fel[0, + oc)(^fe) 



. fc=l 



fe=0 



e"^- exp - ^ (1 - e-"^-) E[e^^^'=l[o,+oo)(^fe)] ■ 



fc=i 



/c=0 



(7.1) 
(7.2) 
(7.3) 



Proof 

Formula ()7.ip is stated in |21| without proof. So, we produce a proof below which is rather similar to one 
lying in |21j related to the maximum functional of the X^s. 

• Step 1. Set, for any (xi, . . . , x„) G M" and a G 6„ (©„ being the set of the permutations of 1, 2, ... , n), 

n / k \ 



U{xi,...,Xn) =^l[0,oo) E 



fc=l 
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and 

y(cr;a:i,...,a;„) =^#Cfc(cr)l[o^oo)( X! ) • 

In the definition of above, tfie permutation a is decomposed into n„ cycies: a ~ (ci{a)){c2{cr)) . . . {cn^ ("■))• 
In view of Theorem 2.3 in j21) . we have the equality between the two following sets: 

{U{a{xi), . . .,a{xn)),cr G ©„} ^ {V{a;xi, . . .,Xn),<y G 6„}. 

We then deduce, for any bounded Borel functions and F, 



E[0(X„)^^([/(ei,...,e„))] = ^ J2 

77 I Z ^ 



In particular, for (j){x) = e*''^ and F{x) = e (where ^ G M and > are fixed). 



E 



,v-y„-,.i:lLi 1[o,+oo)(e ■=! 



fc=ijecfc((T) fc=i Viect((T) // 

^ E exp^M 5] 6-^(#Cfe(a))l[o,oo)( E ^0) 

■cr6S„fc=l L V jGCfcCcr) \jeck(<j) )). 



(Tge„ fc=i 



exp Mm E - ^(#CfeM)l[o,oo) E 



i=i 



Denote by ri{a) the number of cycles of length i in a for any £ G {1, . . . , n}. We have ri{a) + 2r2(cr) 
• • • + nrn{a) — n. Then, 



E[e-— 1 = ^ E n (K^ 



JMX,-1/O[o,oo)(-Yi!) 



E 



N, 



A;i,...,fc^>0: 
fci+2/c2H h'''iA;„— n 



n 



where Nki^...,kn is the number of the permutations ct of n objects satisfying ri(cr) — ki, . . . , r„(cr) = 
this number is equal to 

^fci,....fe„ 



(fci!l'=i)(fc2!2'=2)...(fc„!7i'='.)' 



Then, 



" 1 



fei,...,fc„>0: e=i 

ki +2A;2H \-nkn—n 



E 



,i/iXf-iy£l[o,oo)(Xf) 



• Step 2. Therefore, the identity between the generating functions follows: for \z\ < 1, 



Ee[' 

71=0 



E n 



n>0,A;i,...,fcTi>0: £^1 
fci+2fc2H \-nkn—n 



E 



E n 

fcl,fc2,--->0 £ = 1 



E 



JtiXe-uilio,^l{X() 



n 



k^l 



[]expfE[. 



^iA^^f-Zy^ 1(0,00) (X^) 



JflXi-ul l[0,oo)(^«) 
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expi > JtL|e 

\ri=l 



Z 

n 



The proof of (|7.ip is finished 



• Step 3. 

Using the elementary identity e"^'^^^^ — 1 = (e° — l)lyi(x) and noticing that Tfc = Tk-i + l[o,+oo)(-'^fe): 
we get for any fc > 1, 



gJ'l[o. + cx=)(^fc) _ I 



e" "1 



Now, since Xk — X^-i +S.k where Xk-i and are independent and have the same distribution as ^i, 
we have, for A; > 1, 



E(e*^«i)E(e 



Therefore, 



fe=l k=l 

\ fe=i 



oo 



k=l 



= -^f(.E(e-«0-l)f:E[e 

\ k=0 

By putting (|7.ip into (|7.4p . we extract 

oo 

^E[e'^^'=-^^''l[o,+^)(X,)] z'^ = [e" - (l - zE(e*^«^)) 5(a.,z.,z)] 



(7.4) 



(7.5) 



k=0 



where we set 



/ oo 

S{fj,,iy,z) = expl ^e|^( 



'iM^fc-l'fel[0, + oc)(^fc) 



Next, using the elementary identity 1 ~ ( — exp[log(l — C)] = cxp[— X^fcLi C*^/^] valid for |^| < 1, 



l-zE(e*^«i) =exp( [E{e'^'^')]''^ ^ \ =exp( 



k=l 



fc=l 



and then 



/ oo 

(l-zE(e*^«i))S'(/i,z/,z) =exp ^e[: 



, fe=l 



exp - £ (1 - e-'"^) E[e^'^^'=l[o,+oo)(^fc)] ^ 



fc=i 



(7.6) 



Hence, by putting (|7.6p into (|7.5p . formula (|7.2p entails. 

By subtracting (17.51) from (|7.ip . we obtain the intermediate representation 

oo 

^E[e»^^^-^^'=]I(_^,o)(X,)] = [(e''-zE(e^^«^))5(/i,i.,z)-e'']. 

/£ = 



By writing, as previously. 



- zE{e"'^') =e''cxp -^E(e 



fc=i 
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we find 



\ k=l 
/ oo 

= e''exp( ^(l-e-'''=)E[e^'^^'=l(_,,,o)(^fc)] ^ 
fc=i 



k 



Finally, (|7.3p ensues. 



Lemma 7.2. The following identities hold: 



Proof 

We label the set K as {1, 2, 3, ... , By ()2.5p . we know that the Bfc's solve a Vandermonde system. 

Then, by Cramer's formulas, we can write them as fractions of some determinants: — Vk/V where 



V 



1 

01 



1 



and Vk = 



1 1 1 



q#A'-l 



6* 



1 



By expanding the determinant Vk with respect to its fc*^ column and next factorizing it suitably, we easily 
see that 



Vk = i-ir 



k+l 



92 



7fc+l 



(-1) 



With this at hands, we have 



1 

. . . Ok-i 
■ ■ ■ 01-1 



92 



1 

92 



9f^-^ 



n#K-2 n#K-2 
^k-l ''fc+1 



k£K 



1 

92 



1 

92 



,#K-2 



1 

92 



3#X-2 n*K-2 



'k-1 



'k+l 



1 

92 



q*K-2 



We can observe that the sum lying on the above right-hand side is nothing but the expansion of the 
determinant V with respect to its last row multiplied by the sign (— 1)"*^^"^. This immediately ensues 
that /3#K = (-l)^'^~'nfeeK^fe- Similarly, 



-1 - E Bker^' = 2^ E(-i)'^'^r' 



k£K 



k£K 



^*K-2 



1 

%-l 
92 



1 

dk+1 
Ok+1 



n#K-2 n#K-2 



1 

92 



q#K-2 
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The above sum is the expansion with respect to its last row, multiphed by the sign (— l)'^'^ ^, of the 
determinant V' defined as 

1 ... 1 



V' = 



q#K-2 



Let Rq, Ri, R2, . ■ . , R#K-i denote the rows of V'. We perform the substitution R^k-i ^ 

R#K~i + '^fJ^{~^Yo-iRifK-e where the cr^'s are defined by (|2.3p . This substitution does not affect the 
value of V and it transforms, e.g., the first term of the last row into 

1=2 

Recall that ag — 'Ylii<ki< - <kt<if^K ^^i • ■ • ■ We decompose cr^, by isolating the terms involving 0i, into 



2<k2<---<ki<#K 2<ki<k2<---<ki<#K 



where we set a'^j^ = and a'^ = J22<ki<k2<---<ki,<#K ^kiOk2 ■■ - Okf Therefore, we have 

#K #K #K 

1=2 i=1 1=2 



\keK ) 



The foregoing manipulation works similarly for each term of the last row of V . So, we deduce that 
y = (Efcex^fc) V and finally /3#k+i = {UkeK 0k){EkeK (^f^) ■ ■ 



Lemma 7.3. For a > 0, the Mittag-Leffier functions Ei^a md i?i^ct_|_i admit the following integral 
representations: 



( 1 



r(a) 
1 



l + a;'-"e^ / u"-^e-"du if x > 0, 



1- kl^-^e^ / ' 'u"-ie"du I ifx<0, 



Ei^a + l{x) 



r(a) x^ Jo 

. r(a)|xh Jo 
Proof 

Using the series expansion of i?i,Q, we obtain 



u"-ie-"du ifx>0, 

\x\ 

M"-ie"dM ifx<0. 



E 



(r + a — \)x^ 



■'r(r + a) ^ T{r + a) ''^T{r + a) 

1 ^ ^ r— 1 ^ ^ r— 1 ^ ^ 

1 



^ r(r + a- 1) 

r— 1 ^ ^ 



T{a) 



(7.7) 



(7.8) 



(x + 1 — a) i?i,Q(x) + 



a- 1 



Hence, i?i,Q solves the differential equation xE'J ^ (a;) = (x + 1 — a)-Ei,a(a^)+ p^Q.) - In view of this equation, 
we know that Ei^a{x) is of the form 

i^i,a(a;) = \{x) = \{x) \x\^~" 
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where the unknown function A solves 



i {a)x 



a — 1 
a-l 



if a; > 0, 



This imphes that, for a certain xq > and Aq G M (Aq could be different for x > and x < 0), we have, 
for X > 0, 



A(x) 

r(a) jxo 
and, for x < 0, 

a-l r\-\ 



r(a) 



1 



r(a) 



+ / M"-ie-"du +Ao, 



\x\ 



u°'~^e" du I + Ao. 



Then 



El,a{x) = < 



r(a 

1 



' + f Ao - ^">°" ) x^-°e- ^^" "^ 



r(«) 



1 — Q„a; /-a; 



r(a) 



u"-ie""dw ifx>0. 



+ Ao- 



^ |x|i-"e" 



(7.9) 



t"-ie"du ifx<0. 



[Via) V r(a) ; ' ' r(a) 

Because of the singularity of the differential equation at zero, the initial value at zero does not determine 
xo and Aq. Nevertheless, we know that Ei^a is at zero. So, we need to compute E[ ^(x) for x 7^ 0: 



Ao 



1 



r(a) r(a) 
1 

T{a) fH,,„ 



Xo 

\x\ 



dw (xi-"e^ + (1 - a)x-"e^) + 



r(a) 



du) (|x|i-"e^ - (1 - a)|xr"e^) + 



Via) 



if X > 0, 



if X < 0. 



In order to have a -function at 0, we must have 



Ao 



Xr^ 6 



r(a) r(a) Jo 



e "du and Ao 



x"-ie^o 



r(a) r(a) 



u"-^e"dM. (7.10) 



Putting ()7.10p into ()7.9p yields ()7.7p . Next, formula (|7.8p can be deduced from (|7.7p by simply observing 
that, e.g., for x > 0, 



T{a + 1) 



x"e-^ + / u"e-"du = 



r(a + i) V 



Lemma 7.4. For a G (0, 2), the function Ei ^ admits the following representation. For x > 0, 

sin(a7r) p ^i"" 



El^a{x) = X 



Jo C + 1 



e-"« de 



Proof 



Writing = /J^ e (4+^)" du, we obtain 



00 ^1 — a 



e-^«de-p// e^""de / e~(«+i)"dM = r(2-a)£ 
s + 1 Jo 



r(2-a)^,^^_,^_ 
1 — a 



r(l - a) (^x"-ie-^ - r(a) + u^-^e"" du^ 



u"-ie-" du 
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From this, it entails that 



u^-'e~- du = / |1- e--« + r(a) - x^-^e-^ 

which, by ()7.7p . proves Lemma [7.41 ■ 

Lemma 7.5. T/ie Mittag-Leffier functions Ei i and Ei i are related to the error function according to 

^ + v^e^Erf(VJ) forx>0, 
^ - VHe^'ErffVlaH ) forx<0, 

/TT 

El i{x) = ^ + xe"^' + Ixje'^'Erfda;!) for x e R. 

/n particular, 

1 2 

£;i = ^ +a;e'^ Erfc(|x|) for a; < 0. 
Proof ^ 

Recall that Erf (a;) = ^ e"*' = ^ /p" ^ du and Erfc(a;) = 1 - Erf(a;). By ^J}, we have, e.g., 
for X > 0, that 



E^ i{x) = —= I 1 + y/^e"" / —^du = —= + Va;e"^Erf(V^). 
' 2 ^TT y I- . Ill I . I'TT 

On the other hand, we have, for a; < 0, 



El 1 (x) = — , , , , = — + —7-„ r — E-, i (x'^) + xe^ 

from which we immediately extract the aforementioned representation. ■ 
Lemma 7.6. For a ^ 0, the following equality holds: 

/ ^— ^d«; = — Erfc(a). 

Proof 

Put, for a^O, 



F(a) — a TT dw = — ^ dw. 



We plainly have 



2 . . , I \ 2 2 

e'" dw 



F'{a) = -2a / e"" ^ dw = -2a 1 

= -2a / e""''"'dw + 2aF(a). 
Jo 



+ 1 



So, the function F solves the differential equation F'{a) — 2aF{a) — y/n with initial value F{0) — ^. We 
deduce that F{a) has the form F{a) — G{a) e"^ where the unknown function G satisfies G'{a) — — y^e^"^ 
and G(0) = f . This implies that G{a) = ^/^ e"*' d^ = f Erfc(a) and then F{a) = f e"'Erfc(a). 
The proof of Lemma FTHl is established. ■ 



Lemma 7.7. The following identity holds: for A > 0, s > and x > Q, 

POO 
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Proof 

2 

Put / = 0.3/2'*'' I ^x{s-ij) '^^) With the change of variables ^ = V \sz and u = z + cr/s, we 
have 



2 Jo V"'('*-'^)/'^ ^ 



00 e^"^^^ 



dz dcr 



As e~3s 



2 7o ^ 



3/2 



: du dcT. 



With the new change of variable u = s/ a, it becomes 



/ = 



00 ^-Xs 



: dw du = 



) u—Xsv 



2 Ji Ji vmT^^ 



■ du dv 



Vx 



■ du du 



2 Ji Ji \/uv — 1 

where we set a = |^ and b = 2As. With the change of variable {v, w) — {v, |u + this gives 



w—00 />^) 



Va 

Vab Jw=s^ Jv=i 
1 



: du dw 



dv 



du; 



I . w — a , w — b 

arcsm — + arcsm — 

x./!i+b\ \/w'^~ab y'w'^-ab 



Let us integrate by parts this last integral. First, we have 



_d_ 

dw 



arcsm ■ 



/a{w — b) 



and then, 



X 



w — a , w — b 

— I arcsm — + arcsm 



y/w'^ — ab J {w'^ — ab)y/2w — a — b 

,{w — b) + Vb{w — a) _^ 



\JuP- — ab 



y/a + Vb 



y/w"^ — ab 
: du; 



00 noo 



+ 



±± J si± (w"^ — ab)\/2w — a — b 



dw 



fcL + \/b g + ti 

■ e 



A (u; + Vab)y/2w — (a + 5) xV2 7o {w + ^^ + V^)^ 



■ dui 



/a + vo g+t. 

■ e 



■ dw. 



w^+{^) 
As a byproduct, in view of Lemma 17.61 we have 



I = V2- 



Vb Tie ^^+^^^^(Va + Vb\_Ti 



X V2(Va + Vb) 
Lemma 1 7. 71 is proved. H 



V2 



-e^^^^Erfc — = + VAs 



2^ 



Lemma 7.8. For any integer m < N — 1 and any x > 0, 

poo ^ 

/ e-^"/j.,„(u;a;)du = A-^e-«^^^. (f^J^ 
Jo 

Proof 

This formula is proved in [T3] for 0<m<A^— 1. To prove that it holds true also for negative m, we 
directly compute the Laplace transform of /j.,„(u; .x). We have 
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Let us integrate the function H : z PTI ^ fixed a and M such that 3fJ(a) > and Tlf > on the 

contour Tr = {pe''^ eC:ip = 0,pe [0, R]} (J {pe''^ e C : if e {0,-^),p = R} (J {pe'^f e C : if = 
p e (0, i?]}. We get, by residues theorem, 



oo AI-l foo M-1 



^^ + A 



dz = 2i7r Residue 
iVi 

For M = N — m and a = 0j e'^~x, this yields 
Jo 

Hence, (j2.16p is vahd for m < iV — 1. ■ 



Lemma 7.9. The following identity holds: for any x € M and < s < t, 



C^sin(a;Oe-*«'(^^ e""' du^ d^ 
2t-s 



8 <2(i_s)3/2 



1 - 



Erff-i 



2\l t{t- s) 



Proof 
Set 



A = J fsm{x^)e-'^ e''"' du ) d.^, 

x/tt 2t-s 
B^- — -^xe 



8 i2(<_ 5)3/2 8yit3/2 \^ 2t 

Using the expansion of the sine function, we get 



1-^U--Erff-1 



2 V <(t - s) 



(7.11) 



..|(-ir^f e-d.)de 

In order to integrate by parts this last integral, we search for a primitive of ^2"+3e~*^^. With the change 
of variable C = t£,'^ , we have 



2t"+2 



C+^e-'^dC = - 



2t«+2 ^ A:! 



fc=0 



Then, 



OO 



2^' ^' (2n+l)!i"+2 

n=0 ^ 



4^2 



(n + 1)! 
(2n+ l)!t« 



gr(fc + i) ^ t 



J2n+l 



k=0 



k\ \t 



(7.12) 



On the other hand, by expanding Erf(0 = ^ 2^=0 (-1)"^ = ^ we get 



„ 2t-s 



8 t2(t_s)3/2 
TT 



^ (_l)»^2n (_l)»^2»+2 ^^ 2X 



8x/st3/2 1 ^ n!(4t)" n!(4<) 

V \ 77.^0 n— 



TT ^ \2V t(t - S) 

,1=0 V V V ; 



2n+l ^_-[^yi+lj^2n 

(2n + 1)71! 
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8^2 

COO 
p=0 



2p+l x^P 



(-1)^ 



The computation of the above product of series can be carried out as follows: 



,p=0 



2p+l x^P 



(-1)' 



pi i'it)P jy^^i2q+l)q\yUit-s) 



E 

oo 

E 

oo 

E 

oo 

E 



E 



2p+l 



P,9>0 
p+q=n 



{2q+l)p\q\ \t-s 



E 



2n - 2(7 + 1 / 71 



2(j + l \qj\t-s 



(2n + 2) y 

(2n + 2)y 
^ ^ ^ 2a + 1 



t - s 



(4<)" 
7i!(4t)" 



E 

g=0 



qj\t — s 



n ™2n 



(-l)"a;' 
n!(4i)" 



2q + l \t~ s 



«!(4i)" ^ n!(4(t - s))- 



By inserting (|7.14p into (|7.13p . we derive 



B = 



E 

r!=0 



t- S 



n+1 n 



^ 2q + l 



g=0 



n!(4t)" 



Writing now (^)« = (j^ - 1)« = ELo(-l)'"'^ © ( t^)"'' have 



n (n\ 

2q+l \t- s 



q n 



= E(-i)ME(-ir 



^ q—k 



(DO 



2q+l \t- s 



The sum with respect to q can be calculated as follows: 

n /n— fc\ 



q—k 



q—k 



2q + l \k 



n—k 



9=0 



2q + 2fc + 1 



n — k 

q 



^2q+2k 



(-1)* 



,1 /n 



2 Vfc 



"'0 

= (-1) 



1 /n 



B[ n- k + l,k 



2 \k 

fc ^!r(fc + i) 

2fc!r(n+ §) ■ 



(1-0""'^ 



Plugging (fTTT]) into (fTTe]) . and next (|7T6| into (|7T5l) . we obtain 



B = 



E 



(-1)^ 



8*2 ^/F^^n!(4t)» 

V n— 



E 



< - s 
(-1)" (n + 1)! 



n+l 



(n + 1)! ^r(fc+l) ^ t 



■E 



8*2 f-;; n!(4i)" r(n + 



oo 



4*2 



E(-ir 



n=0 



(2n+ 1)!*" 



fc=0 
n+l 

E 

k=0 



t - S 



r(fc + i) / t 



k\ 



t - s 



„2n 



„2n 
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where we used in the last equahty r(ri + |) = ■ view of (|7.12l) and (|7.18p . we see that both 

members of (|7.1ip are equah A = B. ■ 

Lemma 7.10. We have, for x < 0, 
yijj — x) e Mt-B) fly 

= -2 ^ „ xe +2y^ ^— 1 - TT Erfc --J— zx 







and, /or x > 0, 



?/(?/ + xje dy 







^2^.e--+2V^^^ 1-- e--Erfe 



Proof 

We only produce the proof for x < since the case x > is quite similar. We have 




^ e ''■it j y{y — x)e 4^.(t-a)(!^ t) 



e-tr / (y+ f^'jfy _ lL_£l^ ) e 4=(t-.,w-dy 



e / 2/ — -X \e 4=(t-=)^ dyH xe « / ye iHt-s)y ^y 

J-sx/t \ t J i J-sx/t 



— p it 



P , s(t-s) . .2 t 2 s(t-s)(2s-t) ^ 
4=(t-=)y !___^a;2g-— / g 4.(t-s)!/ dy + 2-i ^xe 

sx/t t J-sx/t 



Integrating by parts, we observe that 



°° . o o2 



y2g 4.(t-.)i^" dy = 2 " xe 4tu-s) + 2 — ^ / e 4»('-=)^"dy, 

sx/t ^ ^ J — sx/t 



we finally get the result. 
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